Driven quantum coarsening 
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We study the driven dynamics of quantum coarsening. We analyze models of M-component rotors 
coupled to two electronic reservoirs at different chemical potentials that generate a current threading 
through the system. In the large Af limit we derive the dynamical phase diagram as a function of 
temperature, strength of quantum fluctuations, voltage and coupling to the leads. We show that 
the slow relaxation in the ordering phase is universal. On large time and length scales the dynamics 
are analogous to stochastic classical ones, even for the quantum system driven out of equilibrium at 
zero temperature. We argue that our results apply to generic driven quantum coarsening. 



Phase transitions are central to condensed matter and 
statistical physics. Initially, emphasis was put on clas- 
sical and quantum equilibrium phase transitions. Later, 
attention moved to non- equilibrium phase transitions in 
which quantum fluctuations can be neglected. These are 
realized when a system is forced in a non equilibrium 
steady state (by a shear rate, an external current flow- 
ing through it, etc.) [l|, Q or when it just fails to re- 
lax (e.g. after a quench) and displays aging phenom- 
ena U 01 • The study of steady states in small quan- 
tum systems driven out of equilibrium 5] has been re- 
cently boosted by their relevance for nano-devices. In 
contrast, the effect of a drive on a macroscopic system 
close to a quantum phase transition is a rather unex- 
plored subject. Some works have focused on non-linear 
transport properties close to an (equilibrium) quantum 
phase transition 0, 0, [§| • Others have studied how the 
critical properties are affected by a drive [ol, [l^, . How- 
ever, a global understanding of phase transitions in the 
parameter space T (temperature), V (drive), T (strength 
of quantum fluctuations), and the role played by the en- 
vironment, is still lacking. Furthermore, experiments in 
2d electronic systems [ij, [13] show interesting features 
in the relaxation toward the quantum non-equilibrium 
steady state (QNESS) but these have not been addressed 
theoretically yet (except for [l^). 

A number of intriguing questions arise in the context 
of driven quantum phase transitions, some of which are: 
How long does it take to reach the QNESS after one of 
the parameters T, V, F is changed? Do the systems al- 
ways relax to the QNESS or, as for classical systems, 
do quenches deep in the T, V, F phase diagram lead to 
aging phenomena and glassy dynamics? What are the 
properties of the latter 'doubly non-equilibrium' dynam- 
ics? Are quantum quenches, obtained by changing V and 
F at T = 0, different from their classical counterpart? 

The aim of this work is to answer these questions 
for a class of analytically tractable models, systems of 
M— component quantum rotors that encompass an infi- 
nite range spin-glass and its 3d pure counterpart mod- 
eling coarsening phenomena. Models of quantum ro- 
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FIG. 1: Non-equilibrium phase diagram for reservoirs with a 
much larger bandwidth than all other energy scales {hujF ^ 
J). Close to the critical point Vc = Vc{T = r_= 0) the 
critical lines are non-analytical (NA). Close to Tc and Fc 
they are power laws with exponents given next to the figure. 
The arrow above Fc indicates that the critical surface is pulled 
up by increasing the coupling to the leads. 



tors are non-trivial but still relatively simple and pro- 
vide a coarse-grained description of physical systems such 
as Bose-Hubbard models and double layer antiferromag- 
nets [lil . The out of equilibrium drive is provided by two 
external electron reservoirs that induce a current flowing 
through the system. In the simplest setting Q each ro- 
tor is coupled to two independent reservoirs. Using the 
Schwinger-Keldysh formalism we analyze the out of equi- 
librium dynamics in the large M limit. We find a phase 
transition, see Fig. 1, between a QNESS {V ^ 0) and 
an ordered phase, we study its critical properties and we 
discuss the effect of the environment. 

The model we focus on is an infinite-range quantum 
disordered system made of N M-component rotors in- 
teracting via random Gaussian distributed couplings, Jij , 
with zero mean and variance J'^ /N. Its Hamiltonian is 

F ^ 

^s-^El^-E-^^^"^"^' n^^MVz. (1) 



Kj 



are the M components of the i-th rotor and L; 



with = ~ihd/dnf, are the M{M - 1) 
components of the i-th generalized angu lar momentum 
operator with Lf = J2i_i<u(^i'^)^ [H, T controls the 
strength of quantum fluctuations; as F ^ the model 
approaches the classical M-component Heisenberg fully- 
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connected spin-glass. In the large M limit it is equivalent 
to the quantum fully-connected p = 2 spherical spin- 
glass [13, llMl ■ A mapping to ferromagnetic coarsening in 
the 3d 0(A/') model can be established in the classical 
and large Af limits 0] . As we shall show, this mapping 
holds for the quantum model as well. Thus, it allows us 
to extend our results also to the 3d ferromagnetic model 
in the large M limit. 

The system is coupled to two independent and non- 
interacting 'left' (L) and 'right' (i?) electronic reservoirs 
in equilibrium at different chemical potentials /xl = /i 
and fiii ^ fi + eV and the same temperature T. R 
and L reservoirs act as source and drain, respectively. 
The details of their Hamiltonian are not important since 
in the small spin-bath coupling we concentrate on only 
the electronic Green functions matter. We focus on 
free electrons with the same symmetric density of states, 
PL(e) = Pb.{^) = centered at e = 0, and with typical 
variation scale around hujp . 

Each rotor is coupled non-linearly to its (double) elec- 
tron bath. For example, for M = 3 rotors we take 
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where 



^aU'^Y ^^11< 
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R] 



HLi 2 '^I'iRl 

iiR(L)i^ '^iiR(L)i are the fermionic operators of 
the R{L) reservoirs, TujJui is the electron-bath coupling. 



chosen to be constant: Tuo 



1. 



hiOc- cr" are the Pauli 
Ai is the fermion label 



matrices (a = 1, M). 7 
inside the reservoirs. 

System and reservoirs are uncoupled at time t < and 
evolve with H = Hs + Hb + Hsb at t > 0. The density 
matrix &i t = g ~ gg ® qj^ ® gji, provides the initial 
condition. ^5, gj^ and gj^ correspond to equilibrium of 
the system at temperature Tq 3> 1, and the L and R 
reservoirs at temperature T and chemical potentials fi 
and /I + eV, respectively. For simplicity, gs is taken to 
be the identity; this choice is equivalent to any other one 
uncorrelated with disorder [l^ . 

We analyze the t > dynamics by using the Schwinger- 
Keldysh formalism yielding a functional-inte gral repre- 
sentation of the Heisenberg evolution 1^ 2^, 21 1. Each 
field carries a ± index associated to the forward and back- 
ward evolution. The action corresponding to Eq. ([1]) is 



a=± 



dt 



2T 



i<j 



The path- integral runs over paths such that nf^{t) — 
M, Viat. One may lift this constraint by using the inte- 
gral representation of the Dirac delta. This amounts to 
introducing auxiliary imaginary fields Xiait) and adding 

= Ea=± f / ^^ait){nlit) ~ M) to the action. 

After expanding the system-leads interaction up to sec- 
ond order in g = LOc|oJF^ integrating out the fermionic 
fields, and taking the large M limit we obtain a 
(Feynman-Vernon-like) action for the rotors. The de- 
tailed computation 2l| confirms that several system- 



reservoir coupling that preserve the 0{M) symmetry and 
the addition of different LL and RR couplings do not 
modify our results qualitatively. In short we obtain 

ab=± i 

Sf,(i, t') = -zabhLol [Gf,(t, t')Gt{f, t) + L^R] . 

The electronic Green functions are Gab{t,t') = 
-i{Tiljait)ipl{t')) with tpa{t),ipi{t) the fermionic fields 
and T the time-ordering operator on the closed con- 
tour. It is convenient to change basis and use retarded, 
G'^ = {Gf_~Gf^)/h, advanced, = (G^+-Gf 
and Kcldysh, G| = i(Gf+ + G^?„)/2, Green func- 
tions. The E's transform in a similar way. For identi- 
cal reservoirs at temperature T and chemical potential /i 
{V — 0), the self-energy components verify the usual fluc- 
tuation dissipation relation of a standard bosonic bath 
= ?icoth(/3?icj/2)3Sf (tj), with 13 = 1/T and 

kg = 1. 

Collecting all contributions the total action, S — Ss + 
S\ + SsB, is 0{MN). Given that the zero-source gener- 
ating functional equals one, one can simply compute its 
average over quenched randomness [1£] and use a saddle- 
point evaluation of the resulting path-integral that is ex- 
act in the large M and N limits. The value of \ia{t) 
at the saddle point is a spatially homogenous function 
\{t). Its time-dependence is determined by the condi- 



tion {nf^(t)) = M with the average taken over S [21 1 



The macroscopic dynamic order parameters are the 
symmetric two-time correlation and instantaneous lin- 
ear response that in the operator formalism are defined 
as MC{t,t^) = ({n,(t),n,(t„)}/2) and MR{t,t^) = 
6{n,{t))/Sh,{t^)\h=Q = -l/h{[n,{t),n,{t^)])9it - t^). 
The field couples linearly to the i-th rotor and the 
last identity is the Kubo formula valid in linear response. 
The exact Schwinger-Dyson equations then read: 



Vit)Rit, t^) = 6{t -~t^)+ j dt" J:R{t, t")R{t", t^) , 

v{t)C{t,t^)^ f dt" ^R{tX)c{t" ,t^) 

Jo 

' ^ dt" ^K{t,t")R{t^,t") , (2) 



with V{t) = n^T-^dl 



-r A(i:), the retarded and 
Keldysh self-energies given by Tiji — -|- J^R and 

-la2 

ft 



Y.K = Sf-HJ^G, and \(t) = -r?V-^dl^C(t,t^ 
t-)+ll dt" [SflXt, t")G{t, t")+j:K{t, t")Rit, t")]. 

In the QNESS, the dynamics are stationary but they 
do not satisfy the fluctuation-dissipation theorem (FDT) 
when V ^ 0. The Lagrange multiplier approaches a 
constant, X(t) — > Aoo, and the linear response satisfies 
a closed equation that once Fourier transformed reads 
R{uj) = -[h^T-^uj^ - X00 + + J^R{to)]-\ The 
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physical solution to this quadratic equation is the one 
that satisfies R{u! —^ oo) = 0. The correlation is given 
by C{uj) = (w)|i?(aj)p and the spherical constraint, 
C{t,t) = 1, implies 



yB 



2ti SEf (w) 



= 1/2 . 



(3) 



The phase transition occurs when R{lo = 0) = J dtR{t) 
ceases to be real, indicating that the stationary condition 
necessary to Fourier transform is no longer valid. Con- 
comitantly, the derivatives of R{uj) in w = diverge and 
hence the real-time response function shows a power law 
decay. This happens when Aoo = Ac = 2 J + (cj = 0). 
Inserting Ac in Eq. we then obtain the equation for 
the critical manifold in the T, V, T space (for a given g 
and huop)- We shall derive the critical manifold for dif- 
ferent reservoirs in full detail in (2l| ; we summarize here 
some of the salient features. 

We first consider g — > after the long-time limit such 
that the asymptotic regime has been established and we 
take hojp much larger than any other energy scale. For 
r = F = we recover the classical critical temperature, 
Tc = J At F = T = we obtain fc = (37r/4)V, 

as for^_the p = 2 quantum spherical model in equilib- 
rium ^16'] and its dynamics coupled to an equilibrium 
oscillator bath [l^. Finally, the critical point Vc on the 
r = T = line is determined by 



de 



2J 



= PL (/i) PR (p.) 



that can be solved numerically and also analytically in 
some special cases. In the large bandwidth limit hutp ^ 
J and we find eVc = 2J for the pl — PR = P symmetric 
case, with = and differentiable at the origin. The 
form of the critical lines are shown in Fig. 1. 

As for finite Tilof we find that Vc varies (contrary 
to Tc and Fc) upon decreasing Tuof/J: the critical line 
Vc{T^T = 0) is re-entrant and, for a single band, eV is 
bounded when the R reservoir is filled. 

When the coupling to the electronic reservoirs, g, is 
finite the critical line in the F = plane remains unal- 
tered but the critical surface on the F direction is pulled 
'upwards' enlarging the low temperature phase for in- 
creasing values of g. This is similar to what was found 
for quantum oscillator Ohmic baths and is due to a spin- 
localization-like effect [3, ■ 

We now turn to the dynamics. Our numerical and ana- 
lytical analysis of Eqs. ([2]) show that after a quench in the 
low-T, weak-F and weak-l^ phase, the dynamics do not 
reach a QNESS 2l|. There is a separation of two-time 
scales typical of aging phenomena First, a stationary 
regime for short time differences t — with respect to 
the waiting-time after the quench, , in which the sym- 
metric correlation approaches a plateau asymptotically 
in the time-difference. Later, an aging regime in which 
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FIG. 2: Dynamics in the driven coarsening regime: numerical 
solution to Eqs. ((U after a quench to T = 0.2 J, V = 0.2 J, 
F = 1, (? = 1 with hujF ~ lOJ. (a) The symmetric correla- 
tion C{t,tw)- (b) The integrated linear response, x(^i^™) = 
J^* dt'R{t, t'} against C, for t^n = 1024 and using f as a pa- 
rameter. The curved part corresponds to the stationary and 
oscillatory regime with {t — t^) /t^ while the straight line 
is for times in the monotonia aging decay of C. 



C depends on the two times explicitly. This behavior 
is shown in Fig. HJa). The plateau value (/ea, so-called 
Edwards- Anderson parameter, measures the fraction of 
frozen rotor fluctuations on timescales much smaller than 
t.u,. The stationary decay depends on all control param- 
eters. (7ea approaches one a.iT^T = V = Q and zero 
on the critical manifold as in a second order transition. 
In the aging regime, the correlation normalized by ^ea is 
identical to the classical one 



-3/2 



(4) 



for < t/tu, < 1. Wc shall prove this result and unveil 
the connection with coarsening anticipated previously by 
exploiting the quadradic form of the full action in the n 
fields. Under the Keldysh rotation (n+,n^) (in,n) 
with m = i(n+ — n^)/h and n = (n+ + n')/2 the 
action is identical to the Martin-Siggia-Rose one for 
a classical Langevin process in a harmonic potential 
~ Kt)6ij)ni ■ Uj. The noise statistics is, how- 
ever, peculiar: because of the quantum origin of the en- 
vironment it has memory, depends on h and satisfies the 
quantum FDT in the V = case. The Langevin equa- 
tions are rendered independent - apart from a residual 
coupling through the Lagrange multiplier - by a rotation 
onto the basis that diagonalizes the interaction matrix 



J, 



eigenvector associated to the eigenvalue J^. The analy- 



and m,. 



with vl the 



sis then follows the same route as in [13, see [23. One 
finds quite naturally that the long-time dynamics cor- 
respond to a Bose-Einstein-like condensation process of 
the M A^-dimensional 'vectors' on the direction of the 
edge eigenvector. The relaxation is then controlled by 
the decay of p(J^) close to its edge. For Gaussian i.i.d. 
couplings p{Jfj,) oc (2J — J^)^/^. This coincides with the 
distribution of the modulus of the Laplacian eigenvalues, 
^j^2p(d/2-i) ^ ^f^2i^2xi/2 d = 3. For this reason all 

models with a square root singularity of the distribution 
of 'masses' J^, as the ferromagnetic rotor model in c? = 3 
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and the completely connected spin glass rotor model, are 
characterized by the same long-time dynamics. Now let 
us show that the aging dynamics are indeed equivalent 
to their classical counterpart. In the ordered phase, tak- 
ing the long tw and t — limits with t/tw fixed (low 
frequency aging regime) the second-time derivatives in 
the efi^ective Langevin equations can be neglected. Fur- 
thermore, only the low-frequency (w <^ ^/ PK) behavior 
of the kernels plays a role in this regime. In this — > 
limit, Yj^{uj) ^ ct G 3?, and S_r(w) is linear. Therefore 
the noise kernels approach a classical Ohmic white-noise 
limit with 'temperature' 



T* = lim Ef (cj)/[29„3Sf (cj) 



(5) 



At y = one gets T* = T. Instead, at T = and 
eV hujF, one has T* — eV/2: the voltage plays the 
role of a bath temperature. This fact has already been 
reported and it is at the root of the derivation of the 
stochastic Gilbert equation for a spin under bias [1^. 
Having argued that the long-time dynamics is governed 
by a classical Langevin equation at temperature T* , it 
is justified that the correlation scales as the one in [17| 
for t/t^ — 0(1), see Eq. a result with two inter- 
esting consequences. In the case of (large M) quantum 
3d coarsening the classical-quantum mapping extends to 
space-time correlations [ilj and proves the existence of 

1 /2 

a growing coherence length ^{tw) oc tj over which the 
rotors are oriented in the same direction and provides a 
real-space interpretation of aging. Moreover, in the same 
long-time regime, the linear response also scales as in the 
classical problem. Therefore, the quantum fluctuation- 
dissipation relation between integrated linear-response, 
X, and symmetric correlation approaches the classical 
one, X ^ ct -|- ((Tea — C)/Tcg, with an infinite effec- 
tive temperature [25|, Tcs oo, as shown in Fig. [H^b) 
(see also [l^, 24 1). In short, the asymptotic aging dy- 
namics are universal, in the sense that the scaling func- 
tions do not depend on T,T,V in the coarsening phase 
and, hence, are equivalent to the classical un-driven ones 
(T = V = 0) [1^. This result mirrors the one obtained 
in [§| for steady state dynamics. 

The environment plays a dual role: its quantum char- 
acter basically determines the phase diagram but the 
coarsening process at long times and large length-scales 
only 'feels' a classical white bath at temperature T* . 
The two-time dependent decoherence phenomenon (ab- 
sence of oscillations, validity of a classical FDT when 
t/tjju = 0(1), etc.) is intimately related to the develop- 
ment of a non-zero (actually infinite) effective temper- 
ature, Tcff, of the system as defined from the deviation 
from the (quantum) FDT [g^. Toff should be distin- 
guished from T* as the former is generated not only by 
the environment but by the system interactions as well 
(Teff > even at T* = 0, llil). Moreover, we found 



netic coarsening (T = 'fixed-point' scenario): after a 
suitable normalization of the observables that takes into 
account all microscopic fluctuations {e.g. Qea) the scal- 
ing functions are independent of all parameters including 
V and r. Although we proved this result through a map- 
ping to a Langevin equation that applies to quadratic 
models only, we expect it to hold in all instances with 
the same type of ordered phase, say ferromagnetic, and 
a long-time aging dynamics dominated by the slow mo- 
tion of large domains. Thus, a large class of coarsening 
systems (classical, quantum, pure and disordered) should 
be characterized by the same scaling functions. It could 
be worth studying carefully systems evolving by barrier 
crossing, a rapid process in which not only the low fre- 
quency behavior of the bath may be relevant. 

We thank C. Chamon, L. Chaput, A. Millis and A. 
Mitra for useful discussions. LFC is a member of lUF. 
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